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We derive a closed form expression for the sum of all the infrared divergent contributions to the 
free-energy of a gas of gravitons. An important ingredient of our calculation is the use of a gauge 
fixing procedure such that the graviton propagator becomes both traceless and transverse. This 
has been shown to be possible, in a previous work, using a general gauge fixing procedure, in the 
context of the lowest order expansion of the Einstein-Hilbert action, describing non-interacting spin 
two fields. In order to encompass the problems involving thermal loops, such as the resummation 
of the free-energy, in the present work, we have extended this procedure to the situations when the 
interactions are taken into account. 

PACS numbers: 



I. INTRODUCTION 

The usual perturbative expansion of the free-energy 
of a gas of massless bosons at temperature T contains 
infrared divergences when three or more loops are taken 
into account. When all the infrared divergent diagrams 
are summed, the resulting expression is finite but non- 
analytic in the coupling constant. This is well known 
in the context of scalar or spin-one gauge fields [J, H, 
3]. The same situation is expected for spin- two gauge 
fields when we apply the finite temperature field theory 
techniques to the graviton gauge field described by the 
weak field expansion of the Einstein-Hilbert action 
One of the main purposes of the present paper is to obtain 
the explicit full result for the summation of these so called 
ring diagrams in the case of gravity. 

In figure [1] we show a typical infrared divergent ring 
diagram containing two insertions of the graviton self- 
energy. It illustrates the two important ingredients in 
the analysis of any of the higher order ring contributions. 
First we need the full tensor structure of the dominant 
high temperature contribution to the graviton self-energy 
(the blob in figure [1]) in the static limit. This is a well 
known quantity which has been studied and shown to be 
gauge independent [5|, |6| . Secondly, we need the three- 
level graviton propagator (the curly line in figure [1} con- 
necting the two self-energies in the ring diagram. Once 
we know these quantities any ring diagram can be ob- 
tained by multiple insertions of the self-energy in a closed 
loop of gravitons, in a rather straightforward manner. 




FIG. 1: The lowest order infrared divergent contribution to 
the free-energy. The curly line represents the graviton propa- 
gator aj 3 and the blob represents the graviton self-energy. 



At this point the tensor properties of the free-propagator 
are essential in order to obtain a closed form and simple 
result for the resummed free-energy. As we will see, if 
the propagator satisfies the following traceless-transverse 
(TT) conditions 

vTD^iJt) = (1.1a) 

k»DH Xa {k) = 0, (1.1b) 

then the sum of the ring diagrams acquire a simple form 
in terms of the three TT projections of the static high- 
temperature limit of the graviton self-energy. 

Taking into account the gauge independence of the 
leading contributions of the ring diagrams, one can 
choose the most convenient gauge for the graviton prop- 
agator. It has been shown in Ref. [7fl that it is possible 
to choose a gauge such that the graviton propagator be- 
comes TT. In principle we could just assume that the 
known gauge invariant result for the graviton self-energy 
could be used safely in the ring diagrams. However, since 
the derivation of the TT propagator is only possible if we 
consider some modifications in the usual Faddeev-Popov 
procedure, which necessarily leads to new ghosts and in- 
teractions, we have also performed the explicit calcula- 
tion of the leading static graviton self-energy and verified 
that the result is indeed the correct gauge invariant one. 
Considering that this calculation involves some rather 
nontrivial cancellations of diagrams, it also constitutes 
an important test of the gauge fixing procedure intro- 
duced in 0. 

This paper is organized as follows. In section II we will 
review the generalization of the Faddeev-Popov formal- 
ism which leads to a traceless-transverse graviton propa- 
gator. We extend the procedure presented in reference 0| 
by taking into consideration interactions of gravitons and 
ghosts more explicitly. In section III we present the calcu- 
lation of the high temperature static limit of the graviton 
self-energy. Having verified that our gauge fixing proce- 
dure yields the correct gauge independent result for the 
self-energy, we proceed, in section IV, to the calculation 
of the sum of all the infrared divergent contributions to 
the free-energy. Finally, in section V we discuss the main 
results. 
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II. GENERAL GAUGE FIXING AND 
FEYNMAN RULES 

In this section we will follow the basic idea of section II 
of Ref. [3] in order to derive the TT propagator. In addi- 
tion to that derivation, we will also obtain the interaction 
vertices involving three types of ghosts and gravitons (see 
also the reference [8j for an analogous derivation in the 
context of spin-one fields were it has been shown that a 
global gauge invariance, analogous to the BRST invari- 
ance is present in the effective action). This involves a 
generalization of the well known Faddeev-Popov proce- 
dure [§, [l(| to the cases when the gauge fixing condition 
is non-quadratic. 

In order to illustrate the main features of the general- 
ization of the Faddeev-Popov procedure, it is convenient 
to consider the following integral over the components of 
a n-dimensional vector h 



where 



Z = J dhexpS{h), 



S{h) = -h T Mh + Si(h). 



(2.1) 



(2.2) 



Later we will associate h with the graviton field h^ v {x); 
the first and second terms Eq. (|2.2[) will be identified 
respectively with the quadratic and interaction terms 
which arise from the weak field expansion of the Einstein- 
Hilbert action. 

Let us now consider the interesting case when S(h) 
is invariant under an infinitesimal transformation of the 
form 



h^h + A(h)9, 



where the operator A(h) is of first order in the derivative 
operator d as well as in h . This symmetry makes the 
integration in Eq. (|2.ip undefined so that we have to em- 
ploy the Faddeev-Popov procedure which leads to the in- 
troduction of a "gauge fixing term" , yielding a quadratic 
term of the form hMh, such that = M 1 is well 

defined. In the context of gauge field theories, is 
the free propagator which will be dependent on the spe- 
cific choice of gauge fixing. In the case of gravity, some 
rather general gau ge fi xing conditions have been inves- 
tigated previously However, in the Ref. it has 
been shown that it is not possible to obtain a graviton 
propagator satisfying the TT conditions given in (jl.ip . 
using the standard gauge fixing procedures. In what fol- 
lows we will present the main steps which are involved in 
the generalized gauge fixing and them apply the results 
to the case of the Einstein-Hilbert action. 

First, we introduce the following factors of "1" in the 
integrand of Eq. (|2.ip 



1= / d6i8(F{h + aA6 x )) det(aFA) (2.4a) 
1= [ d9 2 S(G(h + a A9 2 ))dct{aGA) (2.4b) 



(mr) 



and 

J dpdgexp f-^-fNcfj detN. (2.4c) 

where the N is a "Nielsen-Kallosh" factor [H, [H, [H| , F 
and G are two independent operators, which we will as- 
sume that are of first order in <9, and 5 is the Dirac delta 
function. The use of two operators makes the gauge fix- 
ing prescription more general; the usual Faddeev-Popov 
procedure would introduce only two factors of "1" , which 
corresponds to make the special identification F = G. 

The next step consists in using the infinitesimal gauge 
transformation 



h-> h- a.AQ\. 



(2.5) 



After integrating out the p and q variables and setting, 
for simplicity, N equal to the identity operator, we obtain 

Z = (a) 2 " / d(?i j 6.0 J dhdet(FA)det{GA) 
1 



x exp 



-h T ( M + -F T G ) h 



+Si(h) - h T F T GA9 



(2.6) 



where 9 = 62 — 0\ . 

If we want to avoid the use of mixed propagators which 
would be generated by the last term in (|2.6p , then we also 
have to perform the shift 



M ■ 



1 



E T Q) {F t GA)9, (2.7) 



(2.3) so that (|2.6j) yields 



Z = J d9 J dhdet(FA)det(GA)J 



x exp 



M + —F G ) h 



sAh+\{M + l -eG 



~F (A T G T F) (m + ^F T G 



{F T QA) • 



(f t ga) 0} 



(2.8) 



where we have dropped the infinite normalization fac- 
tors as well as the integration over the gauge orbit J d9\ . 
Since the transformation (|2.7p is not a gauge transfor- 
mation, its Jacobian may not be equal to one. For this 
reason, we have also introduced the Jacobian factor J in 
the integrand of (|2.8p . 

The determinants in Eq. (|2.8p can be exponentiated 
using the standard Berezin integral 



(2.9) 



deti?= / dcdc exp(— c Be) 



3 



where c, c are Grassmann vectors; the first two determi- written as 
nants in Eq. (|2.6p lead to "Faddeev-Popov" like ghosts 
and the field 9 is a "Bosonic" ghost. There would be also 
an extra ghost field associated with the determinant J. 
However, one can argue that in the applications to be 
considered in the present work, these extra ghost fields 
will not contribute. The final expression for Z can be where 



Z = dddcdc dddd dh J exp [S, 



effl 



(2.10) 



J 



S eff ~ 



-h T (m + ^E T 0) h + ^ (h + i (m + ^F t g) (F t GA) 9 



a 



(A T G T F) (m + ±F T G ) (F T GA) 9- c T FAc — d T GAd. 



(2.11) 



r 



Let us now consider the specific example when the gen- 
eral "action" in (|2.11[) describes the behaviour of spin- two 
fields. In this case, the classical dynamics is obtained 
from the Einstein- Hilbcrt action 



Seh — T 



d d a 



-gR, 



where k 2 — 32irG (G has mass dimension 2 - 
Ricci scalar, g = det and 



(2.12) 
d) R is the 



(2.13) 



is the definition of the metric in terms of the graviton 
field hfj, u , which is to be associated with h. It is straight- 
forward to obtain the corresponding expression for the 
quadratic operator M as well as the first interaction 
terms (there will be an infinity number of graviton self- 
interactions) when (|2.12p is expanded in powers of the k. 
From the lowest order quadratic contribution one obtains 
the following tensor components of M 



(¥) 



If 



--{ v ^d a d 13 + ri al3 d^d v ) (2.14) 



can now be identified with the tensor components of A 
in (5m 

Let us now introduce the tensor components of the 
gauge fixing conditions F and G. Following the analysis 
performed in Ref. 0] we choose 



{Q) X,>1V = 



gi rf LV d A +rf LX d v (2.17a) 
g 2 rf lv d x +rf lX 3 v , (2.17b) 



where g\ and gi are two independent gauge parameters. 
It has been shown that this choice is sufficiently general in 
order to make the propagator TT when the limit a — > is 
taken. In addition, the conditions defined by (|2.17p also 
interpolate continuously between other usual gauges such 
as the de Donder gauge in which case g\ = g 2 = —1/2 
and a = 1. 

From the explicit expressions for M, A, F and G given 
respectively by Eqs. (|2~1"4|) , (|2~T6| and ([2~T7|) . together 
with Eq. (|2.1ip and the relevant expressions for Si n t, the 
momentum space Feynman rules can now be obtained. 
Notice also that the multiplication rules for the operators 
F, G and A are such that 



(F T Q) 



(GA)\ = G S ^A^ 



/if, A j 



(2.18) 



(2.19) 



The Einstein-Hilbert action (|2.12|) is invariant un- 
der the space-time dependent coordinate transformation 
x^ — > x^ — 9^ (x) , where 0^ (x) are the generators of the 
transformation. This induces the following gauge trans- 
formation of the graviton field 



nh^ -> Khfiu + {A) llVv x(h)9 > 



(2.15) 



where 



+ »(A^ + W^(«jiM) (2.16) 



{F T GAf\ = F s G S ' a(3 A a 



0,X 



(2.20) 



with analogous relations for the other terms in Eq. 
(|2.11j) . Also, the individual terms in the action are such 
that, for instance, 

h T Adh = J d d xh ll „(M) fnya ' 3 h al3 (2.21) 



and 



c(FA)c = J d d xc x {FA)\c s . (2.22) 
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The momentum space Feynman rules can now be de- 
rived in the usual fashion using the above identifications 
for each term in the action (|2.11|) . (This tedious but 
straightforward task has been done with the help of the 
computer algebra program HIP [15] .) The expressions for 
the graviton and the 9 ghosts propagators are explicitly 
given in Ref. For completeness (and also to define 
the present normalizations and notations) let us briefly 
rederive those expressions. Here we will also give the 
explicit expressions for the c and d ghost propagators. 

The free graviton propagator in momentum space, 
Du V . a p(k), can be readily obtained from the the first term 
in Eq. (|2. 1 1 [) as the solution of 



-F G 



a/3; A7 



(2.23) 



where the subscript "symm" indicates that we have taken 
into account the Bosonic symmetry of the graviton field 
h^y. Also, it is implicit that we have already made the 
Fourier transformations to the momentum space. From 
the symmetry properties under interchange of tensor in- 
dices, it follows that we can parametrize the solution of 
Eq. (|2.23[) in terms of five quantities Cj, % = 1, • • • 5 as 



D 



liv;a.{3 



(fc) 



1 

2fc2 



C'X^ a/sCO, 



(2.24) 



where 



1 



T, 



(k) 



k 2 
1 

- p 



k^kijkoL kp . 



(2.25) 



In terms of this tensor basis, the solutions for C l are 
given by Eqs. (48) of Ref. [7]. They depend on gi, g 2 , 
a and the space-time dimension d, in such a way that 
the TT property (11. 1|) is fulfilled when the limit a — > is 
taken. In this case, the propagator can be expressed as 



D 



TT 

fif,a{3 



1 



PayPal 



1 

d- 1 



(2.26) 



Let us now consider the quadratic ghost sectors of the 
action. From the last three terms in Eq. (|2.11[) one 
can obtain the quadratic and interaction terms for the 
three ghost fields. The quadratic terms can be readily ob- 
tained considering the contribution of the first two terms 
in ()2.16p . From these quantities the ghost propagators 
associated to the fields c M and d^ are given respectively 
by 



and 



DUk) = 



DL{k) 



(2 5 i 



(2 5 2 



2(51 + l)fc 2 



iX^-v) 



2(52 + l)fc 2 



(2.28a) 



(2.28b) 



Similarly the 9 sector of (12. lip yields the following ex- 
pression for the propagator of the 9 ghost 



t IV 



2 

ak 4 



l(d-l) 



1 



4(5i 



1 



8a(d-2) V.91 + 1 52 + l 



kfj, k^ 



(2.29) 



where 



The interaction vertices can also be derived directly 
from Eq. (|2 . . Let us recall that the quantity Si rep- 
resents all the interaction terms, starting with the three- 
graviton vertex, which arise from the expansion of the 
Einstein-Hilbcrt action in powers of n. Some of the ex- 
pressions for the graviton self-interaction vertices have 
been derived before up to the five-graviton vertex [h|. 
Since the argument of Sj in Eq. (|2.1ip has been shifted 
by a 9 dependent quantity, there will also be additional 
interaction terms between the 9 field and the graviton. In 
the figure [3 we show some of the new vertices involving 
the this type of 0-graviton interactions. The numbers 
inside the blobs are meant to indicate when the corre- 
sponding vertex comes from the cubic or quartic terms 
of Si. Here we are only considering the vertices which 
will contribute to the one-loop graviton self-energy. The 
third term in Eq. (|2.1ip also yields new 0-graviton inter- 
actions. In this case, there is only two diagrams which 
are shown in the figure[3J Finally, in the figure[¥]we show 
the cubic and quartic graviton self-interactions, as well as 
the two diagrams involving the interaction of the gravi- 
ton with the two types of Fermionic ghosts. In the next 
section all the vertices shown in the above figures will be 
employed in order to obtain the known gauge invariant 
result for the leading high temperature limit of the static 
graviton self-energy. 



p k^k v 



(2.27) 



(notice that the transversality k^P^ = 0, idempotency 
PjiaPv — Puv an d the trace r/^P^ = d — 1 guarantee 
the TT condition (JTTTJ) ) . 



III. THE STATIC SELF-ENERGY AT FINITE 
TEMPERATURE 

The static graviton self-energy at finite temperature is 
well known quantity and one of the simplest examples 




(a) (b) 




(c) (d) 



FIG. 2: Interactions between the graviton and the 6 ghost 
arising from the second term in Eq. (|2.11|l . The graviton and 
the ghost fields are respectively depicted by curly and wavy 
lines. Graphs (a), (b) and (c) arise from the shifted cubic 
term and graph (d) from the shifted quartic term. 




(a) (b) 



FIG. 3: Interactions between the graviton and the 6 ghost 
arising from the third term in Eq. (|2,11[) . 



exhibiting gauge independence Therefore, it can be 

used as a rather non-trivial test of the gauge fixing pro- 
cedure presented in the previous section. In fact, there is 
an even simpler example, namely the one-graviton func- 
tion. This has been considered in Ref. Q, and used as 
a test of the vertices involving only the interactions of 
three particles. In the case of the self-energy, the contri- 
bution of all the vertices shown in the previous section 
will be taken into account as it can be seen in the figures 
[5] and [5] Therefore, one of the main results of this sec- 
tion will be the verification of the gauge independence of 
the static graviton self-energy, which will be employed in 
the next section, together with the TT graviton propa- 
gator, in order to derive the ressumation of the infrared 
divergences of the free-energy. 

Because of the algebraic complexity involved in the cal- 
culation of some of these diagrams, we have considered 
the special case when the external momentum vanishes. 
Nonetheless, for the dominant high temperature contri- 
bution, this special choice has the interesting property 
of being equivalent to the static limit, as we have found 
explicitly in the early stages of the present investigation. 
More recently, this equivalence of the static and the zero 
four-momenta limits has been proved to be true for all 
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(a) (b) 



(c) (d) 



FIG. 4: Figures (a) and (b) represent the graviton self- 
interactions from the Einstein-Hilbert action and figures (c) 
and (d) represent the ghost-graviton interactions from the last 
two terms in Eq. (|2.11|l . 




(a) (b) 
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-P 



(e) (/) 

FIG. 5: Diagrams which contribute the static limit of the 
graviton self-energy. The curly and wave lines represent 
respectively gravitons and the 6 ghost. The dashed and 
dot-dashed lines represent the two types Fermionic ghosts. 
Graphs (a), (b), (e) and (f) have a symmetry factor 1/2. A 
factor of (—1) is associated with the Fermionic ghost loops in 
figures (c) and (d). 

the thermal Green's functions [l7| (this is not so, how- 
ever, for the long wavelenght limit [181]). 

All the zero momentum diagrams (a) to (j), as well as 
the sum of the diagrams (k) and (1), in figures [5] and [5] 
are such that their integrands have a tensor symmetry 



6 



Tnnnnrv/Tnnririr 




prescription such as 



(ff) 



(h) 




a p 




-P 




OO / OO \ 

J n=-oo J \ n=0 / 



= / (P- 1 p(l + 2((0)) =0. (3.2) 

However, for our present purpose, it would be interesting 
if the quantities 



c 1 = E C I 



(3.3) 



i=i 



happen to be gauge independent even before the sum 
and integration is performed. In order to investigate this 
possibility, let us first write down the results which one 
would obtain when the static self-energy is directly com- 
puted in the deDonder gauge. The calculation is much 
simpler in this case; only the diagrams (a) and (b) as well 
as one of the Fermionic ghost loops contribute to C ! . A 
straightforward calculation yields 



FIG. 6: Contributions involving of the 9 ghosts (wavy lines) to 
the static limit of the graviton self-energy. Except for graph 
(j), all graphs have a symmetry factor 1/2. 



j = d(d-5) _ 1 2 = 4 

16 2' 



0. 



d-3 



C° = -2C d 



(3.4) 



which allows one to parametrize then in terms of the 
basis given in Eqs. (|2.25|) . Therefore, we can express the 
static thermal self-energy as 




(3.1) 



71 — — OO \i — 1 



where I = 1, • ■ • ,10 and 7=11 labels the contributions 
of the diagrams from (a) to (j) and the one from the sum 
of diagrams (k) and (1), respectively. We are using the 
imaginary time formalism [J 0, H| , so that the integrand 
depends on n through the Matsubara frequencies po — 
2imT. 

Once we compute the contributions to integrands 
of II^ Q/3 , the expressions for C\ can be obtained in 
a straightforward manner contracting the integrand of 
II' a/3 with the five tensors T % a ^ and solving the sys- 
tem of five equations. Notice that the momentum inde- 
pendence of the tensors T^ va p and T^ a/3 may require a 



where the second term in C 1 is the only non-vanishing 
contribution which comes from the ghost loop. 

Let us now consider the individual contributions of the 
diagrams in figures [5] and [6] which arise in the case of a 
general gauge fixing. The simplest diagrams are the ones 
shown in figures [5] (c) and (d) . Using the results of the 
previous section, we find that each of these ghost loops 
yield an identical contribution such that 

Cl = C\ = ~, Ci = Cl = 0, i = 2,---,5. (3.5) 

The 9 ghost diagrams (e) and (f) in figure [5] have prop- 
agators and vertices which depend on the three gauge 
parameters as described in the previous section. Despite 
this, the calculation yields a result such that all gauge 
parameter dependence cancels out and we are left with 

Cl+Cl = ^ Cl+Ci = 0, z = 2,-..,5. (3.6) 



Notice that the integrand of one of the ghost loops in Eq. 
(|3.5p cancels with the corresponding #-loop in Eq. (|3.6p 
in such a way that C\ + C\ + C£ + C<1 = -1/2, which 
is the same as the result that one would obtain in the 
deDonder gauge for the single Fermionic ghost loop. 
Let us now introduce the quantities 



C} — C l , i — 1, • • • , 5, 



(3.7) 



7=1 



7 



where C % are given by (|3.4[) . If the gauge independence 
manifests at the integrand level, then 



A' 



11 

E 

7=7 



(3.8) 



should be verified. In the appendix we display the results 
for A 1 and C\ (I = 7, ■ ■ ■ , 11) from these results it can 
be verified that (|3.8[) is indeed satisfied. The expressions 
in the appendix shows how individual diagrams can have 
a very involved dependence on the three gauge parame- 
ters a, g\ and gi. This constitutes a rather non-trivial 
example of the consistence of the general gauge fixing 
procedure presented in the previous section. 

Finally, let us substitute the gauge invariant results 
given in Eq. (|3.4p into Eq. (|3.1|) and perform the sum 
and integration. This straightforward and standard cal- 
culation yields the following temperature dependent ex- 
pression for the static graviton self-energy in d space-time 
dimensions 



i=l 



(3.9) 



where u is the heat bath four velocity and 



C 1 = 
C 2 = 

c 3 = 
c 4 = 



djd - 3) 2< 4 ~ d ) %2 

(d-i)r(^i) 
-c\ 



dC\ 



-Gp d , 



(3.10) 



C 5 = -d(d + 2)C 1 ; p d = 



_ r(d)C(«0w 



T a , 



with r and £ respectively the the Euler's and Riemann's 
functions. From these expressions we obtain, for d = 4, 
the results given in Eqs. (3.8) of Ref. [f| (notice that in 
the present work we have labeled the tensors T l in a dif- 
ferent order, so that the constants Cj of Ref. [|| are such 
that ci = C x /k 2 = p 4 /12, c 2 = C 3 = 0, c 3 = C 5 /k 2 = 
-2p 4 , c 4 = C 2 /k 2 = -pi/12 and c 5 = C 4 /k 2 = p 4 /3 
when d = 4). It is also possible to include contributions 
of other thermal particles, such as scalars or fermions. 
Such contributions will only modify all the C 1 by a com- 
mon integer factor which counts the number of degrees 
freedom associated with each field. 



IV. RESUMMATION OF INFRARED 
DIVERGENCES IN THE FREE-ENERGY 

The sum of all the infrared divergent contributions to 
the free-energy can be represented graphically as 




(4.1) 



where the quantity II stands for the static graviton self- 
energy computed in the previous section. It is implicit 
in (|4.ip that we are considering only the zero mode con- 
tribution of the graviton propagator (denoted by curly 

lines), so that k 2 — — \k\ 2 . The reason for this is because 
each photon propagator in (|4.ip introduces a factor 



1 

k 2 



1 



(27m T) 2 



(4.2) 



so that the zero mode (n = 0) of each individual ring di- 
agram in (|4.ip yields an infrared divergent contribution, 
when d = 4, from the momentum integration J d 3 k ■ ■ ■ . 
In scalar theories, as well as in gauge theories of spin one 
fields, it is well known that the sum of all such infrared 
divergent contributions yields a finite result which is non- 
analytic in the coupling constant ■ In what follows we 
will employ the results of the last two sections in order to 
investigate how a similar result may be achieved in the 
case of spin two fields. The details of this analysis show 
explicitly that the use of the TT graviton propagator al- 
lows one to obtain an explicit and compact result. 

Since we are going to employ the TT graviton prop- 
agator, the computation of the right hand side of Eq. 
(|4.ip becomes much simpler if we express the static self- 
energy in terms of its traceless and transverse compo- 
nents. However, at finite temperature there are three 
distinct TT tensors tensors which may depend on the 
momentum k and the heat bath four velocity u. In the 
static case, when k ■ u = fco = (we are adopting the rest 
frame of the heat bath), these tensors can be written in 
d dimensions as 



1 



d-2 



(4.3a) 



l 



d- 1 

T^ Xa (u,k)-T^ Xa (u,k) (4.3b) 
kn.kjj 



(d - Vju^Uy + - JL — 
(d - l)u\u a + 



k 2 
k\k a 



k 2 



VXa 



where 



VflV ~ V^v 



- u,,u v - 



k^k v 



,(4.3c) 



(4.4) 



For d = 4 the above expressions reduces to the static 
limit of the corresponding ones given in Ref. 0. No- 
tice that the sum of the three traceless transverse ten- 
sors coincides with the numerator of the propagator in 
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Eq. (|2.26[) , which represents the only traceless transverse 
tensor available at zero temperature, so that we can write 



jyTT _ _j_ (rpA , rpB , rpG \ 



(4.5) 



In terms of its TT components the dominant contribu- 
tion to the static self-energy can then be expressed as 



Similarly, in the case of the higher order graphs we obtain 

I /d(d-3) {CA) n + {d _ 2){CB) n + ((?C) r 



H*l a 



(4.13) 



n 



aT™/3 — C / T^ ly Q/3 + C 4 T 4 V a p + ■ ■ ■ , (4.6) 

I=A,B,C 



where the ellipsis represents terms which are orthogonal 
to the TT tensors and the constants C 1 and C 4 can be 
solved in terms the ones given in Eq. (]3 . 10[) as follows 



C A = C B = 2C 1 



C 4 = C 4 



d-1 
C 5 
d-1 



(4.7a) 
(4.7b) 

(4.7c) 



Notice that the dependence of each individual TT ten- 
sors T 1 on the momentum variable k = (0, k) is canceled 
when all the terms in (|4.6[) are taken into account so that 
the result becomes identical to the one in Eq. (|3.9p . This 
apparent unnecessary complication is very convenient in 
order to perform the calculation of the ring diagrams in 
(|4.ip . This happens because the TT tensors in (|4.3[) are 
not only traceless and transverse, but also enjoy of some 
very important properties. First, they are idempotent so 
that 

t Iv\s tI a p = T ^ a 0i [I = A,B,C). (4.8) 
They are also orthogonal 



Their "norm" is given by 

rpA rpAV-v a<5 _ d(d — 3) 
J^^vxs 1 — 2 

Finally, the tensors T A and T B also satisfies 



T A xs u x = T* x5 u x u 5 = 0. 



(4.9) 

(4.10a) 

(4.10b) 
(4.10c) 

(4.11) 



Using these properties, as well as Eqs. (14.5[) and (|4.6p it 
is straightforward to show that the integrand of the first 
ring diagram in Eq. (|4.ip is given by 

r-jTT /ii^/ii v\ yrterm. T)TT AifTi ^2 f2 Tjierm. 
I /d(d_3) (C - A) 2 + {d _ 2)(C - B) 2 + { - Cy 



fc 4 



where we have used the zero mode condition k 2 = — \k\ 2 . 
Notice that the components of the self-energy which are 
not traceless and transverse drops out in the final result. 
This is consistent with the fact that the TT components 
are physical. 

We have now all the ingredients to compute the sum 
of the infrared divergent contributions to the free-energy 
in Eq. (|4.1|) . From Eq. (14. 13)) we can see that the sum of 
the ring diagrams is composed of three similar structures, 
each one having the same form as — ~Y^ =2 {—x) n /n — 
log(l + x) — x, so that the free-energy can be written as 



n(T) 



T 



2 (2tt 



d{d - 3) i(c A ) 



(d-2)I{C B )+I{C c )} 



where 



1(c) 



d d ~ 4 k 



log ( 1 + - 

\ \k\ 2 \k? 



(4.14) 



(4.15) 



is a familiar integral which arises also in the context of 
scalar or vector fields and it can be done in a closed 
form. Performing the d — 2 angular integral and using 
the change of variable z = \k\/ ^/c, we obtain 



dz 



,d-2 



r(¥) jo 

(4.16) 

Using integration by parts, and performing the resulting 
integral yields 



7(c) = - 



4(7TC) £ 



(d-i)r(^i) 



oo z d ~ 4 

dz- 



z z + l 



4r (¥) 



(d-l)(d-3) 



(ttc)" 



(4.17) 



Substituting (|4.17[) into (|4.14[) we obtain the following 
result 



n(T) 



2r(^)T 



(d- l)(d-3)(27r) 



d-1 



d(d-3) (7T - A) ^ 



+ (d-2)(7rC s )^ + (irC c Y 



(4.18) 



(4.12) Finally, using the Eqs. (|3.10|) and (|4.7[) this expression 
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can be written as 



n(T) 



2 r(^) 



(d-l)(d-3){2n) d - 1 

2(5-rf) 7r 5 T t r(d + i)C(d) 
r(ifi) 

d(d-3) fd-3 
d-l 



+{d-2) 



d- 3 



1 



d-l 

(d- 2)d(d + 2)\ fd-3 



2(d-l) 
< (GT d - 2 )^ T d 



d-l 



(4.19) 



DISCUSSION 



In this work we have investigated the possibility of ex- 
tending the general gauge fixing procedure of Ref. Q 
in order to take into account interacting spin two fields. 
As an explicit example, we have computed the dominant 
one-loop contributions to the thermal self-energy of the 
graviton, and verified that it agrees with the known gauge 
invariant result. Then, using a decomposition of the self- 
energy in terms of the three traceless and transverse ten- 
sors which arises at finite temperature, as well as the TT 
graviton propagator, we were able to obtain a closed form 
expression for the sum of the infrared divergent contri- 
butions to the free energy. 

We note that the above gauge invariant result for the 
graviton thermal self-energy satisfies a simple Ward iden- 
tity [l!| (not just a BRST identity) which is a conse- 
quence of the fact that the ghost self-energies are sub- 
leading at high temperature. Consequently, on dimen- 
sional grounds, the thermal self-energies associated with 
the c M and d^ ghost fields must be of order GT d ~ 2 k 2 , 
whereas the 9^ ghost self-energy would be proportional 
to GT d ~ 2 k 4 , where the factor GT d ~ 2 is dimensionless. 
The presence of these powers of k in the thermal ghost 
self-energies ensures that infrared divergences will be ab- 
sent in ring diagrams involving ghost propagators. This 
justifies the neglect of such diagrams in the evaluation of 
the leading infrared divergent contributions to the free 



energy of spin-two fields. (A similar behavior occurs also 
in the case of the free-energy in QCD [l|.) 

The result presented in Eq. (|4.19p has some interest- 
ing features which we would like to stress. First, for odd 
space-time dimensions it is a real and singular function. 
On the other hand, for even space-time dimensions, it is 
a finite and non-analytic function of GT d ~ 2 as one would 
expect for a non-perturbative quantity. However, in this 
case it acquires an imaginary part. For instance, for d = 4 
the third term inside the curly brackets, which can be 
traced back to the T component of the self-energy, is 
equal to (— 7/3) 3 / 2 . As a result, one would conclude that 
the gravitational C-mode is unstable, since the imagi- 
nary part of the free energy is connected with the decay 
rate of the quantum vacuum [20| . However, a detailed 
investigation shows that the graviton self-energy, which 
is proportional to GT 4 , is of the same order as the so- 
lution of the Einstein equation for the curvature tensor, 
when the thermal energy momentum tensor is taken into 
account. Therefore, by consistency, one should also take 
into account the curvature corrections in the analysis of 
instabilities of gravity at finite temperature. These cor- 
rections [1, @ have the effect of adding some extra contri- 
butions to the self-energy in such a way that the C-mode 
contribution to £1{T) would change the third term of the 
curly bracket of Eq. (|^7Tg|) to (-7/3 + 5/27) 3/2 , which is 
still imaginary. This term may be related to an imaginary 
value of a thermal Jeans mass [fl @], which reflects the 
instability of the system due to the universal attractive 
nature of gravity. 
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APPENDIX 

Here we display the expressions for the quantities in- 
troduced in Eq. 



J 



(3d -7) ( 9l 



92) 



8 [(d - l)(<?i - g 2 ) 2 + 2a(d-2)( 9l + l)(.g 2 + 1)] 



(la) 



Ci = 



(lb) 



10 



Gl 



(d 2 - lid +22) (gi- ff2 ) 2 



8 (d-2) [(d-l)( 3l - 52 ) 2 + 2 a (d-2)(. gi + l)(.g 2 + l)] 



G* = 



(1c) 
(Id) 



(d- 3) 2 (.gi - ,g 2 ) 2 



10 



2 (d-2) [(d-l)( 3l -.g 2 ) 2 + 2 a (d-2)(. 9l + l)(. 92 + l)] 



Gl, =0 



(le) 
(If) 



(51 - 52 ) 2 4 (d - 3) ( 5l - , 92 ) 2 + a (d - 1) (d - 2) (<?! + 1) ( ff2 + 1) 



4 [(d - l)(5i - ,g 2 ) 2 + 2a(d - 2)(. 9l + 1)(. 92 + 1)]' 



G 7 2 = 



(2a) 
(2b) 



cl = 



(d 2 - lid +26) (pi - , 92 ) 2 



8 (d - 2) [(d - l)( 5l - , 92 ) 2 + 2a(d - 2)((/i + 1)(. 92 + 1)] 



G Q 2 = 



(d- 5) 2 (. 9 i -5 2 ) 4 



C 



8 [(d - l)(. 9 i - , 92 ) 2 + 2a(d - 2)(</i + l)(fl2 + l)r 



( 5 i - ff2 ) 2 2 (d - 3) 2 (51 - 52 ) 2 + a (d - 2) (d 2 - 7d + 14) (51 + 1) (g 2 + 1) 
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2 (d - 2) [(d - l)(5i - 5 2 ) 2 + 2a(d - 2)( 5i + l)(.g 2 + 1)]' 



G11 = 



(2c) 
(2d) 

(2c) 
(2f) 



A 3 



(d - 3) (51 - 92? 



4 [(d - l)(5i - 52)2 + 2a(d - 2)(5i + l)( ff2 + 1)] 



G 3 = 



(3a) 
(3b) 



Go 



(d-3) (d-4) (5i-. 92 ) 2 



4 (d-2) [(d-l)(5i-5 2 ) 2 + 2a(d-2)(5 1 + l)(5 2 + l)] 



G 3 = 



(3c) 
(3d) 



r 3 - 
°io — 



(d - 3) 2 (51 - g 2 f 



2 {d-2) [{d - l)(5i - 92? + 2a(d - 2)(</i + l)( ff2 + 1)] 



(3e) 
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ct, = o 



(3f) 



(d - 3) ( 5l - g 2 f 2 ( 5l - g 2 f + a(d-2) ( 9l + 1) (g 2 + 1) 



2 [(d - l)(. 9 i - , 92 ) 2 + 2a(d - 2)(</i + 1)(. 92 + l)] z 



(4a) 



C 7 4 = - 



(d-5) (.9i - .92)' 



7 ~ 2 [(d - l)(5i - g 2 ) 2 + 2a(d - 2)(. 9 i + l)(g 2 + 1)] 



(4b) 



ci = 



(d-6) (d-3) (.91-.92)' 



4 (d-2) [(d-l)(.9i-.9 2 ) 2 + 2 a (d-2)( 3 i + l)(.9 2 + l)] 



(4c) 



(d-3) (d-5) ( 9l -g 2 ) 4 



4 [(d - l)( 5 i - , 92 ) 2 + 2a(d - 2) (.91 + l)(</ 2 + 1)Y 



(4d) 



r 4 _ 
<-^m — 



(d - 3) (.91 - ,9 2 ) 2 (d - 3) (,9i - . 92 ) 2 + a (d - 2) (d - 4) (.91 + 1) (52 + 1) 



(d-2) [(d - l)(.9i - . 92 ) 2 + 2a(d - 2)0/1 + 1)0/2 + l)] 2 



(4e) 



°11 



A 5 = 



C 7 5 



(d - 5) (.91 - .92) 



cl 



■y5 
"10 



C 



11 



2[(d- 




-. 92 ) 2 + 2a(d 


- 2)(,/i + 1)0/2 + 1)] 






(rf-3) 2 ( 5 i- 


52) 4 


2[(d- 


1)0/1- 


-</ 2 ) 2 + 2a(d- 


2)0/1 + 1)0/2 + I)] 2 






(2d-6) (</i - 


- 52) 2 


[(d- 


i)0?i- 


</ 2 ) 2 + 2a(d- 


2)0/1 + 1)0/2 + 1)] 






Cf = 








(d-3) 2 (51- 


S2) 4 


2[(d- 


1)0/1- 


- 52 ) 2 + 2a(d- 


2)0/1 + 1)0/2 + I)] 2 






(d - 3) 2 (51 


~52) 4 


[(d- 


-i)0/i- 


- ff2 ) 2 + 2a(d- 


-2)( 5 i + l)( ff2 + l)] 2 






(2 d-6) ( 9l 


-.92) 2 


[(d- 


- l)(ffl 


- ff2 ) 2 + 2a(d- 


-2)(.9i + l)(.92 + l)]' 



(4f) 

(5a) 

(5b) 
(5c) 

(5d) 

(5e) 

(5f) 
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